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Abstract 

We derive reduced finite dimensional ODE models starting from one dimensional 
lubrication equations describing coarsening dynamics of droplets in nanometric poly- 
mer film interacting on a hydrophobically coated solid substrate in the presence of 
large slippage at the liquid/solid interface. In the limiting case of infinite slip length 
corresponding in applications to free films a collision/absorption model then arises 
and is solved explicitly. The exact coarsening law is derived for it analytically and 
confirmed numerically. Existence of a threshold for the decay of initial distributions 
of droplet distances at infinity at which the coarsening rates switch from algebraic 
to exponential ones is shown. 

1 Introduction. 

Dewetting processes of a liquid polymer film of nanometer thickness interacting on a hy- 
drophobically coated solid substrate attracted an intensive research during last several 
decades, see. e.g. a review in jlj. In general, such processes can be divided into three 
stages. During the first stage a liquid polymer film is susceptible to instability due to 
small perturbations of the film profile. Typically such films rupture, thereby initiating 
a complex dewetting process, see e.g. [2|-|4|]. The influence of intermolecular forces play 
an important part in the rupture and subsequent dewetting process, see e.g. [1, @] and 
references therein. Typically the competition between the long-range attractive van der 
Waals and short-range Born repulsive intermolecular forces reduces the unstable film to 
an ultra-thin layer that connects the evolving patterns and is given by the minimum of 
the corresponding intermolecular potential, i.e. the film settles into an energetically more 
favorable state, see 0, 0]. The second stage is associated with the formation of regions of 
this minimal thickness, bounded by moving rims that connect to the undisturbed film, see 
e.g. 0-Ei ■ 

In this study we are interested in the third and the last stage of the dewetting process, 
namely the long-time coarsening process that originates in the breaking up of the evolving 
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patterns into small droplets and is characterized by its subsequent slow-time coarsening 
dynamics, which has been observed and investigated experimentally by Limary and Green 



12l . Il3| . They show that during the coarsening the average size of droplets increases and 
the number of droplets decreases. The coarsening mechanisms that were observed in such 
films are typically subsequent collapses of smaller droplets and collisions of neighboring 
ones. During collapse the size of a droplet shrinks in time and its mass is distributed in the 
ultra-thin layer. In turn, collisions among droplets occur due to the mass transfer through 
the ultra-thin layer between them that causes a translation movement of them, droplet 
migration, eventually leading to the formation of new droplets. A numerical example of 
the coarsening dynamics in two-dimensional films is shown in Fig. [2j 




Figure 1: Plots of intermolecular pressure Il e (/i) (blue) and potential function U £ (h) (green) for 
e = 0.1 



Besides intermolecular forces and surface tension at the free surface of the film the 
dewetting of polymer films on hydrophobic substrates also involves such boundary effect as 
slippage on a solid substrate 14|. Recently in Munch et al. 



15 closed-form one-dimensional 



lubrication equations over a wide range of slip lengths were derived from the underlying 
equations for conservation of mass and momentum, together with boundary conditions for 
the tangential and normal stresses, as well as the kinematic condition at the free boundary, 
impermeability and Navier-slip condition at the liquid-solid interface. Asymptotic argu- 
ments, based on the magnitude of the slip length show that within a lubrication scaling 



there are two distinguished regimes, see [15 



These are the well-known weak-slip model 



d t h = -d x M(h)d x (ad xx h - U e (h)) 



'l.F 



with M(h) := h 3 + bh 2 and b denoting the slip length parameter; and the strong-slip model 



Re (ed t (hu) + d x (hu 2 )) = vd x (hd x u) + hd x (ad xx h - U E (h)) - — (1.2a) 

ed t h = - d x (hu) , (1.2b) 

respectively. Here, u(x,t) and h(x,t) denote the average velocity in the lateral direction 
and the height profile for the free surface, respectively. The positive slip length parameters 
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Figure 2: Numerical solution to ()1.2a|) - (jl.2bj) with e = 0.1, /3 = 2.5, Re = 1 showing an example 
of a coarsening process (collapse of the 4th small droplet and collision of 2nd and 3rd ones) in 
the array of five quasiequilibrium droplets. 



h and (5 are related by orders of magnitude via b ~ f] 2 P, where the (small) parameter rj, 
<t) < 1, refers to the vertical to horizontal scale separation of the thin film. 

The high order of the lubrication equations fll.ll) and (jl.2ap — (11.2bp is a result of the 



contribution from surface tension at the free boundary, reflected by the linearized curvature 
term ad xx h with parameter a > 0. A further contribution to the pressure is denoted by 
H £ (h) and represents that of the intermolecular forces, namely long-range attractive van der 
Waals and short-range Born repulsive intermolecular forces. A commonly used expression 
for it 0, [l8| is given by 

2 2 

n ^) = ~ Ti with < e < 1. (1.3) 

It can be written as a derivative of the potential function U e {h) = U{h/e) (see Fig. [TJ) 
where 

The parameter < e 1 is the global minimum of U e (h) and gives to the leading order 
thickness of the ultra-thin layer. Below we often use notation for the pressure and flux 
functions function 

p(h) := ad xx h - U E (h), j{h) = hu (1.5) 

Re (d t (hu)+d x (hu 2 )) and ud x (hd x u) in fll.2ap - fll.2bl) . with Re, v > denoting the Reynolds 
number and viscosity parameter, represent inertial and Trouton viscosity terms, respec- 
tively. 

Additionally, the weak-slip and the strong-slip models contain as limiting cases three 
further lubrication models. One of them is the no-slip model, which is obtained setting 
b = in the weak-slip model: 

d t h=-d x (h 3 d x (d xx h-ii £ (h))y (i.6) 
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The second one is obtained from the strong-slip model in the limit (3 — » oo and describes 
the dynamics of suspended or falling free films: 



Re (ed t (hu) + d x (hu 2 )) = ud x (hd x u) + hd x (ad xx h — II e (/i)) (1.7a) 
ed t h = -d x (hu), (1.7b) 

For the third limiting case the slip-length parameter /?/ is of order of magnitude lying in 
between those that lead to the weak and the strong-slip model, i.e. b <C f3j <C (3. The 
corresponding intermediate- slip model is given by 

d t h = -d x (h 2 d x (d xx h - U e (h)) ) . (1.8) 

It can be obtained by rescaling time in (11. ip by b and letting b — > oo or by rescaling time 
and the horizontal velocity by in (ll.2al) - (11.2bl) and taking the limit /3 — > 0. Existence of 
weak solutions to (jl.2ap - (11.2bl) and (ll.7al) - (11.7bj) and rigorous conver genc e of the former 



ones to the classical solutions of (II. 8p as f3 — > was shown recently in [171 ] . 



As in 17|] we consider systems (ll.2ap - (ll.2bp and (ll.7al) - (11.7bp on a bounded interval 



(0, L) with the boundary conditions 

u = 0, and d x h = at x — 0, L, (1.9) 

whereas equations (II .11) . (II .61) and (ll.8p with 

d xxx h = 0, and d x h = at x = ±L. (1-10) 

Both (II. 9p and (11.101) incorporate zero flux at the boundary and as a consequence imply 
the conservation of mass law 



1 f L 

— j h(x, t) dx = const, Vt > 0. 
L Jo 



Within the context of thin liquid films one of the first studies of the coarsening dynamics 



can be found in Glasner and Witelski [18, 19 . The authors considered the one- dimensional 



no-slip lubrication model ( II. 6p with ( ll.lOp . They confirmed numerically existence of the two 
coarsening driven mechanisms, namely collision and collapse. One of the typical problems 



considered in [18|, [19J was the calculation of the coarsening rates, i.e. how fast the number 
of droplets decreases due to coarsening in time. Often in order to identify the characteristic 
dependence for coarsening rates one needs to model very large arrays of droplets (around 
10 4 ). But due to the presence of the ultrathin-layer of order e between droplets the problem 
of numerical solution for any lubrication equation becomes very stiff in time and demands 
high space resolution as the number of droplets increases. Therefore, there exists a need for 
further reduction of lubrication models to more simple, possibly finite-dimensional ones. 

Basing on the observation that solutions of lubrication equations describing coarsen- 
ing dynamics stay in time very close to a perturbed finite combination of quasistationary 
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droplets and can be therefore parameterized by a finite number of parameters, namely po- 
sitions and pressures of drops, in [18|, |l9| for the first time a reduced ODE model describing 
evolution of the latter ones on the slow time scale was derived from the lubrication equation 
(II. 6p . Using this reduced model the authors derived also the corresponding coarsening law 
in the form 

n(t)~r 2 /\ (1.11) 

where n(t) denotes the number of droplets remaining at time t. Later, analogous reduced 
ODE models from lubrication equation (11.11) with a general mobility M(h) = h q , q > in 



one and two dimensional case were derived and analyzed in |20l . |2JJ . An step to a rigorous 



justification of these models basing on a center manifold approach was made recently in 22 
For the case M(h) = h the coarsening law (11.111) was justified rigorously in 23j using the 



gradient flow structure of the corresponding lubrication equation. The work of |20l . 121 
concerns migration of droplet. There it was shown that the direction of the migration of 
droplets governed by (II. ip with a general mobility M(h) = h q , q > is opposite to the mass 
flux applied to them. Moreover, for g < 2 the driving coarsening mechanism is collapse 
of droplets that is due to the mass diffusion in the ultra-thin layer between droplets and 
similar to Ostwald ripening in binary alloys, see 24j-|26| . Note, that also in the no-slip case 
q = 3, i.e. one described by (II. 6p . as was shown in [l9[ the coarsening rates even for the 
systems coarsening solely due collisions obey the law (II. lip . 



Recently, in Kitavtsev and Wagner [27( it was shown that the coarsening dynamics 



of quasistationary droplets governed by (ll.2ap - fll.2bp with sufficiently small Re number 
is driven also by collapse and collision. There reduced ODE models analogous to that 
one of 18|, |l9] were derived for system (ll.2ap - fll.2bp and its limiting case ( ll.7ap - fll.7bp 
as well. In contrast, to the case of (II. ip it was found there that the coefficients of the 
strong-slip reduced ODE model depend explicitly on the slip length (3. In particular, there 
exists a critical length (3 cr = 0(e) such that the migration of droplets proceeds in the 
direction of the applied mass flux for (3 > (3 cr and opposite to it for [3 < (3 cr . Moreover, it 
was shown that for moderate and large the driving coarsening mechanism switches from 
collapse to collision of droplets. Basing, on these observations it was conjectured and shown 
numerically in (2?J that the coarsening rates for systems fll.2ap - fll.2bl) and fll.7al) - fll.7bp 
can be remarkably different from ones for (II. lip . 



In this study we continue the research initiated in [27] ■ Our aim here is to derive explicit 



coarsening laws for the dynamics of droplets in the strong-slip and free film regimes, i.e. 
governed by lubrication system fll.2ap - fll.2bp and its limiting case fll.7ap - fll.7bp . The 
missing point in 27] was a derivation of flux representation between interacting droplets 
for moderate and large slip lengths which was important for closure of the derived there 
reduced ODE models. Therefore, inspired with the matched asymptotics technique applied 
in Glasner 20j to the lubrication equation fll.ll) . we present in section 2 a new closed form 
derivation of reduced ODE models for fll.2al) - fll.2bp and fll.7ap - fll.7bp that incorporate the 
explicit flux representation for all < < oo. 

In section 3 we concentrate on the reduced ODE model corresponding to fll.7al) - fll.7bl) . 
i.e. on the regime of free films characterized by the infinite slip length = oo. In this case 
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migration and subsequent collisions of droplets dominate completely collapse component 
of the coarsening dynamics. Therefore, we look only at the migration subsystem of the 
derived reduced ODE model such that droplet pressures are kept constant during evolution 
of droplets and updated only after each subsequent collision event. We observe further that 
for a special initial data this migration subsystem can be solved explicitly while its solution 
represents subsequent collisions of iV — 1 droplets with the largest last one. Therefore, we 
call it as an exactly solvable collision/absorption model. It turns out that the coarsening 
law for this model depends only on the initial distribution of the distances between droplets 
and can be derived analytically. Finally, we derive the continuous counterpart of the 
coarsening law proceeding to the limit N — > oo. 

In section 4 we consider several examples of initial distributions of distances between 
droplets and show that the corresponding coarsening rates depend only on the distribution 
decay at infinity. Moreover, for an explicit family of distributions decaying as 1/ x 1+a with 
a > we show existence of a threshold at a = 1 at which the coarsening rates switch from 
algebraic to exponential ones. 

In section 5.1 we justify the derived hierarchy of the reduced models by numerical 
comparison of their solutions to ones of the initial PDE system f ll.2ap - fll.2bj) and its limiting 
cases fll.Sp and fll.2ap - fll.2bj) . We observe that the deviation between them stays 0(e) 
uniformly in time. Besides we compare solutions of the collision/absorption model from 
section 3 with those of the full reduced ODE system for the case (3 = oo. Finally, in section 
5.2 we check numerically the derived coarsening law for the collision/absorption model in 
the case of finite N and its continuous counterpart. 



2 Derivation of reduced ODE models. 

We consider a solution to fll.2ap - fll.2bp which stays close in time to a union of iV + 1 
droplets, which precise characterization to be described below. Similar to the derivation of 
reduced coarsening models for the classical thin film equation in 20] we distinguish three 
regions in our matched asymptotic analysis. 

• Droplet core (DC) region: This region corresponds to droplets and is composed of 
the union of disjoint intervals (Xi(t) — Ri(t),Xi(t) + Ri(t)), so that Xi{t) and Ri(t) 
are the center and the radius of the 2-th droplet, i = 0, N. The dynamical points 
Xi(t) ± Ri(t) are called contact line points and are defined through the relation 

h(X l ±R l )=eH\ (2.1) 

where H* is the global maximum of U'(H) with function U(H) defined in fll.4p . We 
expand 

Ri = Ri,o + zRi,i + Xi — Xifi + eXi t i + ... (2.2) 

and denote 

• dR ■ dX 

R = — — , X = ——. 

at at 



6 



Contact line ( CL ) region is a microscopic internal layer around the contact line points 
where h and x scale like e. Here we employ the moving rescaled spatial coordinate 

z= flft)-!*-*^ da = i djtt dt = d T -- £ d z (R oT X ), (2.3) 

where in this section the sign =p corresponds to two CL regions around the points 
X =F i?, respectively. Accordingly, by definition ( 12. ip we have h(z = 0) = eH*. 

Precursor layer (PL) region is the complement (— L, L) \ Uj(Xj — Ri,Xi + Rj). In 
this region h scales like e. 



we 



The main goal is to determine the evolution of Ri(t) and Xi(t). To do so as in [20 
propose self-consistent asymptotic expansions in each of three regions and connect them via 
matching conditions. Corrections to the leading order base solutions solve linear equations, 
and Fredholm-type solvability conditions will yield information about the dynamics. 

Let us first consider the motion of the i-th droplet with % G 1 , . . . , iV — 1 . For a time 
being we skip below the subscript i. Let us start with the CL region. Here the solution to 
fll.2ap — dl.2bp is expanded as 

h = eHi + e 2 H 2 + u = eU 1 + e 2 U 2 + ... 

We will also use the induced expansions 

P = P + eP x + J = e 2 J 2 + e 3 J 3 + ... 

for the pressure and flux functions defined in (jl.5p . The corresponding leading order system 
in e in this region is given by 

d z (ad zz H 1 -U'(H 1 )) = 0, 
d z H 1 (R TX ) = Td z (H 1 U 1 ). 

Integrating the last system and using matching conditions to the DC and PL regions 

d z H x ->■ 0, d zz H x ->• 0, Hi ->• 1 as z — >■ — oo 

d zz Hi -> 0, #i — ► +oo as +oo (2.4) 



one obtains 



In particular, 



l^) 2 = ^(^-[/(l), (2.5a) 
Ui = -(I-^-PoTXo)T^ (2.5b) 



lim U x = -(R tX ), lim t/j = T-M-oo). (2.6) 
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Next, in the DC region we expand the solution as 

h = h + sh\ + e 2 h 2 + u = u + eu\ + e 2 u 2 + ... 
and correspondingly pressure as 

p = p + ep! + ... 
In turn, the leading order system in this region is given by 

ahM^M - u j = o, 

-d x (h uo) = 
Integrating the system and using the matching condition 

h (X TR) = (2.7) 

one obtains its solution in the form 

ho = R7m {R2 ~ {x ~ x{t))2) ' MoEE °- (2,8) 

Correspondingly, the leading order pressure is given by 

Po = TT7^- ( 2 - 9 ) 



In the PL region we expand the solution as 

h = eh\ + e 2 h 2 + u = eu\ + e 2 u 2 + ... 
and correspondingly pressure and flux as 

p = Po + ep x + j = e 2 j 2 + 5 3 j 3 + ••• 
The leading order system in this region is given by 

hMu'ih)) = o, 

d t hi = -d x (hiui) 

Integrating the system and using the matching condition h\ (X =p R) — 1 one obtains 

hi = 1, ui = j 2 = const (2-10) 
For the next order corrections h 2 ,u 2 in PL region one has the system 

hid x (U"(hi)h 2 ) = 

d t h 2 = -d x (hiu 2 + h 2 ui) 



From the first equation and (12.101) one obtains 

d xx h 2 = d xx p = and j 2 = ~(3d x p . 



(2-11) 



Proceeding further in the expansion in the CL region for the second order corrections 
H 2 , U 2 one obtains the system 



ud^H&Ut) = H 1 ((Td gg H 2 -U"(H 1 )H 2 ), 

o = rt-ioj-^+M; 



(2.12) 



Let us introduce a linear operator 



" H ' 




Hxd z {ad zz H 2 - 


U"(H X )H*)) 


U 




d z (H 2 (R -X - 


- Ux) - H X XJ 2 



Formal adjoint operator to C is given by 
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V 





-ad zzz (H l9 ) + U"{H 1 )d z (H l g) - d z V(R^ - X - U, 

d z vH l 



The kernel of it contains two linear independent functions 



9i 




" 1 " 




92 




' l/#i " 


. Vl . 







5 


. v<1 . 








(2.13) 



To derive necessary Fredholm-type solvability conditions for the system (12.121) we multiply 
the first equation in (I2.12p by g 2 and integrate it on (— oo, +oo) to obtain 



P (+oo) - Po(-oo) 



dziHxdzUx), 



where we have used that in the CL region 

P = U'\H X )H 2 - ad zz H 2 . 
Substituting in the previous expression (I2.5bj) one obtains 

P (+oo) - Po(-oo) = -v/(^J 2 (-oo) + P =F ^o) 
where a constant integral I is given by 

1 



J_s (Mi 

„ Hi z \ Hx 



dz 



35(3 + v 7 ^)' 



(2.14) 



(2.15) 



(2.16) 



and can be effectively calculated from ( 12. 5 aft (see Appendix A). Formula (I2.15P is an analog 
of Gibbs-Thomson boundary condition and shows that the pressure experiences a jump 
at the CL region. Note, that this is a first considerable difference between the coarsening 
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dynamics driven by (jl.2aj) — (jl.2bp and (II. ip . In contrast to (I2.15P as was shown in 18] the 
pressure is constant through the CL region in the case of (II. ip . 

Next, multiplying the first equation in fl 2 . 1 2 [) by and integrating it on (— oo, +oo) 
one obtains 

+ / H 1 d z (ad zz H 2 -U"(H 1 )H 2 ). 



= vE x d z XJ x 

Integrating further three times by parts and using (I2.5bp . (12.141) one arrives at 

d z H t 







—v- 



-( T J 2 (-oo) + R tX 



-H X P Q 



+oo 



ad zz H x H 2 



+oo 



- ad z H 1 d z H 2 

Using the matching condition (12.41) and additionally 

d z H 2 —> const as z — > — oo, 

d z Hi — > d x h , d z H 2 ~ d x h\ + d xx h z, Hi ~ hi + d x h Q z as +oo (2.17) 
one arrives at 



{HiPo 



a d z Hi ( + oo ) d z H 2 ( + oo ) . 



The last expression again using (12. 4p and (I2.17P implies 



cd xx h Q 



a(d x h d x hi 



-P(+oo), 

P(-oo) - P(+oo)hi(X =f R). 



(2.18) 



Note, that the first relation in (I2.18P is consistent with already derived ( I2.8p - (l2.9p . whereas 
the second one is new. 

Finally, let us consider the system for the first order corrections hi, Ui in the DC region 
which has the form 



dh ■ dh ■ 
r-n-o 7T~ A o 



OR dx 
Let us introduce a linear operator 



= vd x (h d x ui) + ah d xxx hi - ui/(3, 
-d x (h ui) 



(2.19a) 
(2.19b) 



" h ' 




u 





isd x (h d x ui) + ah d xxx hi - ui/f3 
-d x (h ui) 



Formal adjoint operator to £ is given by 



9 




V 





ud x (h d x g) - § + h d x v 
-crd xxx (h g) 
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The kernel of it contains two linear independent functions 



9i 




1 



92 
1'2 



1 

C x dr_ 

JX /3ho 



(2.20) 



To derive necessary Fredholm-type solvability conditions for the system fl2.19ap -f l2.19bj) we 
multiply (I2.19bj) by v±, integrate it and using the matching condition (12. 7p obtain 



R = 0. 



(2.21) 



In turn, multiplying the right hand side of the second equation in (12.121) by v 2 and inte- 
grating it on (X — R, X + R) one obtains 







-X 
X 



X+R 



dh, 



x-r 9x 

X+R 



v 2 dx + 



X+R 



d x (h Ui)v2 dx = h Q uiv 2 



X-R 



X+R 



dh 



X-R ^X 

2X R 1 f x+li 



v 2 dx = X (h v 2 ) 
ui dx. 



x+R 2X R 1 



X-R. 



X-R 
X+R 

U\ dx 

X-R 



X+R 



hoU\d x V2 dx — 



X-R 



X-R 



In the last equality we used that h ~ 0(R — \x — X\) and v 2 ~ log(i? — \x — X\) as 
x — > X =p R. Next, using (I2.19ap and integrating three times by parts one arrives at 



2X R 



X+R 



u\ dx 



X-R 



X+R 



ud x (h d x ui) + ah^d^Jix dx 



X-R 



= [vh d x ui + ahod^hx - adxh^d^ + ad xx h Q hi 
Let us note that from fl2TT9aj) - fl2.19bp and (J22IJ, (EH) it follows that 



X+R 

X-R 



r) h 



Ul = Xq, 
X 



(3ah 



(2.22) 

(2.23a) 
(2.23b) 



Hence, by (|2.23aj) the first term in the square brackets in (12.221) vanishes. By f!2.23bp one 
has d xx h\ ~ log(i? — \x — X\) as x — > X ± R. Due to this and fl2.7p -f l2TS|) the second term 
in the square brackets in (12.221) also vanishes. In turn, due to the matching condition 



h 1 (X= F R) = H 1 (0) 



(2.24) 



and (I2.8P the fourth fourth term in the square brackets in (I2.22p vanishes. Therefore, 
relation (I2.22p reduces to 

2X R 



- \ad x h d x hi] 
11 



X+R 
X-R 



(2.25) 



At this moment let us introduce back the droplet subscript i = 1, N — 1 and denote 
by Jj the flux j 2 in the PL region between % — 1-th and i-th droplets. Combining (12. lip 
with (12.151) and (I2.2ip one obtains that Jj is constant and satisfies 

T Q Pj - Pi-i - 2VJJ - uI(X L0 + X,_ 1>0 ) 

■h = P j • (2.26) 

In the last formula we introduced two more notations: the constant leading order pressure 
inside i-th droplet 

Pi = — \=- (2.27) 
Ri\/3a 

according to (12. 9 p and the distance between the neighboring DC regions 

di = Xi — — Ri — Ri-\. 
From (I2.26P one obtains an explicit expression for Jc 

J i = p Pi ~ Pi - 1 ~ U jfi'° R + jti - lfi \ i-h-,N. (2.28) 
di + Ivlp 

Next, from f|27Z5|> . the matching conditions P^gjl . fF2T2^) and equations(J2IE>D, (ESD one 
obtains the leading order equation for i-th droplet position evolution 

Xlfi = ~2R~TW^ (Jl+1 + Jl) - 



Substituting in the last expression the flux representation (12.281) . denoting 

di 

and using (I2.27P one obtains 



di = A, (2.29) 



v P t ((P l+1 -Pi)-Iu{X l+1 + Xi) , (Pi-P^-Iv&i + X^ 

-X-i,0 — 



2/(V3aP) + HvPi \ d i+1 + 2 di + 2 

for % = 1,...,N- 1. (2.30) 

In turn by (I2.2ip and definition (I2.27P one has 

Pi,o = 0, for i = 0,...,JV. 

The derived ODE system describing the leading order in e evolution of pressures and 
positions of N + 1 droplets will be closed if we additionally prescribe that the first and the 
last droplet do not move, i.e 

Ao,o — ^Xn,o = 0, X = 0, X N = L. (2-31) 
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The condition (I2.3ip corresponds to the situation when one extends the array oniV + 1 
droplets from interval (0,L) to an infinite array on the whole real line R by reflection 
around the points x = and x = L. It stays also in agreement with boundary conditions 
(USD and ffTTOl) 

Let us point out that the evolutions of pressures is slower than one of positions and 
proceeds on the order e. One can potentially obtain it by going further in the expansion 
of the solution to fll.2ap - fll.2bl) . while an easier way is to derive it from the conservation of 



droplet volume as was done in [18|, l20| for the case of equation (II. ip . Namely, the volume 
of the i-th droplet V{ is changing due to the difference of the fluxes in the surrounding it 
PL regions. Using (12.81) and (I2.27P one obtains 

4v4 3 • 2 

Substituting in the last expression the flux representation (12.281) and denoting 
one obtains 

eP 41 = - C P i +1 ~ Pi - -C t ( X l +1 - Xl - Xi J Xi A , i = i, N - 1. 

Iv\d i+1 +2 di + 2 J \ d i+1 + 2 di + 2 J 

Finally, combining the last expression with (I2.30p and (I2.3ip the closed ODE system for 
the leading order evolution of positions and pressures in the array of iV + 1 droplets takes 
the following form: 

^ = P / (Pi+i - Pi)- Mji+i + Xi) | (P. - Pi-i) - HXj + Xj. X) 




2/{V3cP) + 21vP, \ d i+1 + 2 di + 2 

P t = Clf Pi+i-Pi _ P t - Pi-A _ f X 1+1 - Xi - X *J X ^ 1. ;■-■.] V -l:(2. 

Iv \ d i+l + 2 di + 2 ) \ d i+1 + 2 di 

and 

P X + 2C 1 J^— = ~\ X 1 = 0, 

d x + 2 lv d x + 2 

P N -2C N -J^ = ^ P 1~ Pn -\ X n = 0. (2.33) 
d N ^ + 2 Iv d N _ x + 2 

Let us consider certain limiting cases. In the case (3 — > oo the limiting system for evolution 
of pressures and positions has the form 

C 

Pi + CjpQ+i — Xi-x) = — (Pi+i — 2Pj + Pi-i), 

Lv 

X i+1 -2X i + X i .. 1 = Pl+1 ~j Pi ~\ for i = 1,...,N-1; (2.34) 
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and 

A + CaX 2 = ^l(p 2 - Pl )^ x 1 = 0, 

Pn — CnXn-i = — ^—(Pn ~ Pn-i), X n = 0. (2.35) 

Lv 

Next, rescaling the time by f3v and proceeding to the limit /3 — >■ the limiting system for 
evolution of pressures and positions takes the form 

Pi+l ~ Pi Pi ~ Pi-1 



i+1 

di+i d{ 



and 



Pi = Ct 

Xi = + for i = l,...,N-l; (2.36) 



P N = -2C n Pn ~/ N -\ X n = 0. (2.37) 



d 



N 



Note, that the last system coincides with one derived in [20[ for the intermediate-slip 



equation (11.81) in the one-dimensional case. This stays in agreement with the fact that 



(II. 8p is the limiting case of (ll.2ap - (11.2bj) as — > as was shown in 15|, [17j]. Finally, note 



that after time rescaling by (3v taking limits v — > oo or v — > results again in ( 12. 34[) — (I2.35j) 
and (I2.36p - (l2.37p . respectively. This is also natural, because (ll.7al) - (ll.7bp and (11.81) are 
the limiting cases of (jl.2ap - (ll.2bl) as well as v — > oo or v — > 0, respectively. 

Let us summarize the algorithm for simulation of coarsening dynamics in large arrays 
of droplets using the derived reduced ODE models. Starting with an array of N + 1 
droplets after each subsequent coarsening event (i.e a collapse of one droplet or collision 
of two droplets) one can model the coarsening process further by reducing the dimension 
of the model by two and solving the reduced ODE model with the updated initial data. 
Practically, as in [3] we say that a collapse event occurs at a moment when pressure of 
one droplet increases a certain threshold, namely when 

27 

P > 0.5P max (e), with P max {e) := — — . (2.38) 

2bbe 

Then we take the final pressures and positions for remaining droplets from the previous 
run of the reduced ODE model as initial conditions for the next one. In the case of collision 
in Glasner and Witelski [19] was suggested that coarsening event occurs when the distance 
between two colliding i-th and % + 1-th droplets becomes smaller then a certain threshold 
5 = 0(e), i.e. when 

di < 5, (2.39) 
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After the collision we calculate the position and the pressure for the new formed droplet 
by formulas 



Xi^ new — l/2(X i+1 — Ri + \ + Xi — Ri), 

( 1 1 v 1/2 

Pi,new — y~p2 p2 J " (2.40) 

The last formula for Pi, new is based on the observation that mass of the new droplet is to 
the leading order in e given by the sum of the masses of the collided droplets (see flit). 
In section 5.1 we compare solutions of the derived reduced ODE model (I2.32I) - (I2.33I) with 
those of the initial PDE system (ll.2al) - (11.2bl) and show that the former ones provide high 



accuracy 0(e) also after subsequent coarsening events. 

3 An exactly solvable collisions/absorption model 

Let us consider the limiting case of infinite slip length = oo, namely the ODE system 
(I2.34p - fi2.35p describing coarsening in free films. As pressure evolution proceeds on a slower 
time scale then that one of positions as e — > let us consider only migration of droplets. 
Namely, we investigate the zero order system 

X n = X N = P i = 0, for i = 0,...,N, 



X i+1 -2X i + X i _ 1 = Pt+ \j Pl ~\ for i = 2,...,JV-l; (3.1) 

As will be justified numerically in section 5.1 for given e, T > one can choose initial 
data with sufficiently small -Pj(O) <C 1, i = 0, 1,...,N such that the difference between 
solutions to (JSU) and f !2TMj) - fl2351) stays uniformly 0(e) for all times t G (0,T). Note that 
for such initial data there is no other constraint on the location of X(0) other than that 
di(0) should be larger then the collision threshold introduced in f)2.39p . 

Moreover, for certain initial data one can solve (13. ip explicitly. Indeed, if 

Pi(0) = p, for i = 0, 1, N - 1 and P N (0) = p with 1 > p > p (3.2) 

then the solution to (13. ip is given by 

Xi(t) = X0) + -£t, for i = 1, N - 1; X = 0, X N = L, where B = (3.3) 

In this and the next sessions for convenience reasons we redefine the notation for the 
distances between droplets from the previous section as 

di(t) = Xi(t) - X_i(t) for i = 1, N - 1 and d N (t) = L - X N ^(t) - R N (t) - R N -i(t) 
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and call usually di(t) as the distance of the i-th droplet at time t. Using this notation one 
can rewrite the solution (13. 3p in the following form 

(U(t) = di{0) + —t, for i = 1, ...,iV — 1 

d N (t) = d N (0) - B{N ^ 1) t for t e (0,T C ), where T c = ^ (3.4) 

Note, that T c denotes the time proceeding until (N — l)-th droplet collides with the largest 
last one. Iterating ( 13 .4p one observes that first N — 1 droplets collide one after another 
with the last one like rings in the famous rubber band toy. Due to (13. 3p all droplets except 
the first and the last ones move to the right. The last droplet consequently absorbs the 
neighbor droplet, while the distance between them is uniformly distributed between the 
remaining droplets. Therefore, the distances of the remaining droplets at the collision time 
T c are given by 

d i (T c ) = d i (0) + d N (0)/(N-l), i = l,...,jV-l. (3.5) 

Writing the solution to (13.11) in the form (I3.4p is convenient because one can substitute 
di(T c ) as the initial distances for the modelling of the next collision event. 

Note, that due to (13. 5p distance monotonicity is preserved in time for solutions (13.41) . 
i.e. if di(0) > d m (0) for some l,m 6 1, ...,N then di(t) > d m (t) for all times t > 0. This 
property allows us basing only on a given initial distribution of the distances in the array 
of droplets to derive the coarsening laws analytically for solutions to (13.11) considered with 
(13. 2p and additional assumption 

l»p»P N . (3.6) 

This assumption prescribes that the last droplet is much larger then others and allows us 
to simplify further the dynamics by assuming that its pressure P/v remains constant in 
time. In turn, this implies that the coarsening dynamics in this case depends solely on the 
evolution of droplet positions without change of their pressures after subsequent collisions. 

Indeed, let initial distances be prescribed by k G N families such that there are i m 
distances in m-th family (1 < m < k), all of them are equal to d m and 

di > d 2 > ... > 4, k + i2 + ... + ik = N (3.7) 

holds. Additionally, let us allocate these k families in the initial configuration so that the 
distances between droplets non increase coming from the first to the last droplet. Then 
due to the distance monotonicity property this ordering will be preserved in time, i.e. first 
the members of the family k will be absorbed by the last droplet, then those of k — 1-th 
one and etc. Moreover, the distances in each family will stay equal for all t > 0. This 
implies that for initial data satisfying (I3.2p . ( 13. 6 j) and (13. 7p all collision times are uniquely 
determined having given k and the set {d m , i m }, m — 1, k. Therefore, using the explicit 
solution (13.41) holding between subsequent collision events the coarsening law can be derived 
analytically by a recursive procedure. 

Indeed, let us fix an index 1 < m < k and look at the moment when all families with 
the indexes m + 1, k and also I — 1 members of the m-th family have been absorbed for 
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some given 1 < I < i m . Let us calculate the time t(n) needed for the absorption of the l-th 
member with n denoting the remaining number of droplets after the latter event. Using 
(13 .4p one can easily calculate by recursion that 



i-i 



n + l ~ ^ t(n + r)B \ (n + l)d m 



r=l 



where by d m we denote the distance in the m-th family at the time when (m + l)-th family 
has been absorbed. From (13 .Sp one can obtain the total time needed for the m-th family 
to be absorbed T m in the form 

T m = % ( N - ± i,) ± - ' ■ . (3.9) 



r-n,+ i / r=l ^ Y2 p =m+l' l P r 



In turn, using again (13 .4p one recursively finds 

~ dpi p _ Ep=m+i dpip 

«m — at _ sr k ~ ~~ at _ V^fe ~ + r ' 

p=m+l iV l^p'=p % P iV l~ip'=p % V 

Substituting the last expression in (13. 9p one obtains 
T m = — Nd m + ^(dp ~ d m )ip ^2 



J N - V i - r 

p=m / r=l Z-^i p=m+l P 



Therefore, the total time needed for all families up to m-th to be absorbed is given by 

Let us now derive the continuum version for the discrete coarsening law in (I3.10P pro- 
ceeding to the limit iV — > oo and k — > oo. Suppose we are given a probability density 
function f(d) on (0, +oo), i.e. 

r+oo 

/ f(x) dx = 1, f(d) > and /(d) = if d < 0. 

Defining dj = iAd for i G N U {0} and a fixed Ad < 1 we approximate f(d) by a piece- wise 
constant function f a (d) as follows. 

f a (d) = f{d i+ i) for d e [di, d i+1 ) and / (d) = if d < 0. 

Accordingly to this approximation suppose we are given an array of iV + 1 droplets with 
N ^> 1 such that the number of droplets with the distances lying in the interval [d iy d i+ i) 
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is equal to [Nf(d i+ i)Ad]. As before we suppose that droplets are allocated so that the 
distances between them are non increasing and (13. 2p . (13. 6p hold. Then using (13.101) one 
obtains 



m / k \ Nf(d p )Ad 

m / k \ Nf(dp)Ad 



Nl, ^ As 



= EjK + EK-W £ TT^ +0(Ad,l/iV) 

P =0 \ p'=p / s=l/N 1 Z^ p '= p +1 J \ u p)^ a b 

Taking the limit N — > +oo in the last expression and introducing 

TAd{d) = hm 

one obtains 

P =o ^ V P '= P / ^ 1 - £ P '= P +i W) Arf - s 

(3.11) 

Applying the Taylor expansion to the last integral in (13. lip one finds 



f ^ )Ad ds _ /(rfp)Arf 

s=0 i-E^ =p+ i/K)a^-^ ~ l-E^=pW)Ad 



Inserting this into (13. lip one obtains 



B P =o\ P '= P / 1 -J2 P '= P f(d' p )Ad 



Finally, taking the limit Ad — > and introducing 

T(d) = hm 

V 7 Ad^O N 

one arrives at 

m = 5 f (* + f<» - x)m *) i-i/U dI (312) 

Introducing function n(cf) as the relative number of droplets with initial distances larger 
or equal d, i.e. as 

n(d) = 1 - / f(x)dx (3.13) 
18 



one obtains from (I3.12p that 



T(d) = U!{- X + £ {V - X)n ' Mdy )T^fTy 

1 f d 

= — / n(x) In 
Jo 

The last expression provides an exact coarsening law, i.e. it tells what time T(d) is needed 
until all droplets having initially distances smaller then d are absorbed by the last large 
droplet. 

In appendix B, we show that the discrete coarsening law (13.101) can be recovered back 
from (13.141) if the initial distribution f(x) has the form 

k 

f(d) = J2^(d-d m ), (3.15) 

m=l 

i.e. if initial distance distribution is represented by k e N families as in (13. 7p while the 
number of droplets iV — > oo. Moreover, we justify the connection between (13.141) . (I3.10P 
and the starting ODE system (13. ip as well numerically in the section 5.2. 



n[x) 



n(d) 



dx. 



(3.14) 



4 Examples of coarsening rates 

a) We consider an explicit family of initial distributions f(x) and show that depending 
on their decay as x — > +oo the coarsening rates reproduce all possible algebraic decays. 
Moreover, there is a certain threshold after which the decay becomes exponential. Namely, 
let us consider a family 



From (I3.13P it follows that 



x 



l+a 



dx A a 



x 



n[x) 



l+a rgl+a 



with a, A > 0. 



x 



Substituting this in (I3.14p one obtains 



T(d) 
T{d) 



olA 



B(a - 1) \a-l 

2 



1-a 



- 1 



+ a In 



if a 7^ 1, 



A 
B 



In 



/2 + ln 



if a = 1. 



Combining (14.3b|) and (14 ,2p one obtains the exact coarsening law for the case a = 1 



n(t) = exp 1 - VI + 2Bt/A 



(4.1) 
(4.2) 

(4.3a) 
(4.3b) 

(4.4) 
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In the case a^l one obtains from (I4.3al) and (14. 2 p 



T(n) 



aA 



1 



a-l 

n « — 1 



B(a - 1) VI - a 
For the latter exact law one obtains the following asymptotics 



ln(r^ 



n(t) ~ < 



tB(a 



exp 



aA 

tB(a - 1) 
aA 



if a < 1 



if a > 1 



as t — > oo. 



(4.5) 



Therefore, from ()4.4p -( H~5"j) one finds out that for < a < 1 the coarsening rates are 
algebraic at least for large times, while at a = 1 they become exponential and stay so for 
a E (1, +oo). 

b) Consider f(x) = exp(— x). Substituting it in (I3.13P and f)3.14p . consequently, one 
obtains the exact law ^ 

T(n) = —{n — 1 — ln(n)). 
B 

Thus, in this case the following asymptotics holds 

n{t) ~ exp(— Bt) as t — > oo. 



(4.6) 



c) Consider a Gaussian distribution f(x) = exp (— x 2 ). In this case by ( 13 . 1 3[) one 

has n(x) = erfc(x). Substituting it in (I3.14p one obtains 



T(d) 



1 

B 



d r ex 



n(y) dy 



n'(x) 



n\x) 



dx 



B^ 



(1 -exp(-x 2 ))exp(-x 2 ) 2 

— 2x exp(— x J 



J x °°exp(-a; 2 )o?t 



dx 



^(-c + o( e M-J>))+l 



exp (—a; 



J x °° exp (-x 2 ) dt 



dx 



BV^ 



-C - \nn(d)) + 0(exp(-d 2 )), 



where constant C ~ 0.74. Therefore the following asymptotics holds 

n{t) ~ exp{— C — By/nt] as t — > oo. 



(4.7) 



and the coarsening rates show an exponential decay as in the example b). 

d) Finally let us show that the coarsening rates for large times depend only on how fast 
the initial distribution f(x) decays as x — > +oo and not its behavior for moderate x. In this 
and the next example we consider non-monotone distributions having a local maximum 
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2 4 6 8 10 12 14 16 18 20 2 4 6 8 1 1 2 14 16 18 20 

X X 

Figure 3: Initial distributions in the example d) (left) and the example e) with a = 2 
(right) 

at x > 0. Consider f(x) = (1 — x) 2 exp(— x) (see Fig. [3]) with n(x) = (1 + x 2 ) exp(— x), 
correspondingly. By (13. 14ft one obtains 

BT(d) = / (1 + x 2 ) exp(— x) ln(l + x 2 ) dx — / x(l + x 2 ) exp(— x) dx 
Jo Jo 

f d 

- (ln(l + rf 2 )-d) / (1 + x 2 )exp(-x)dx 
Jo 

[1 + x 2 ) exp(-x) ln(l + x 2 ) dx - 7 + 3(d - ln(l + d 2 )) + 0(exp(-d)).(4.8) 
The first integral in the last expression can be estimated as follows. 

fd pd 

j (l + x 2 )exp(-x)ln(l + x 2 )(ix < ln(l + d 2 ) j (1 + x 2 ) exp(-x) dx 
Jo Jo 

= ln(l + rf 2 )(3 + 0(exp-d)) 

Combining this with (I4.8P one obtains 

T(d) = ^ + o(d) 

and hence the following asymptotics holds 

9 



o 



n(t) ~ ^1 + -^ 2 j exp(-Bt) as t ->• oo. 

Therefore, the coarsening rates show an exponential decay as in the example b). 
e) Consider distributions 



a 



a 



- [(1 — x) 2 exp(— x) + 1/(1 + a;) 1+a ] with a > 0, a ^ 1. (4.9) 
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They have a local maximum at x > and a decay ~ l/x 1+a as x — > oo (see Fig. [3]). 
Correspondingly, one has 



1 + a 

Substituting it in ( I3.14p one obtains 

" d r rx 1 n'(x) 



[(1 + x)~ a + a exp — + x 2 )] 



BT(d) 



n(y) dy 



(d 



n(x) 

1_a - 1) + a (3 - exp(-d)(3 + d(2 + d))) 



x) 



1 — a 

1 - xf exp(-x) + (1 + xy^ a 



x 



dx. 



(1 + x)~ a + aexp(-x)(l + x 2 ) 
The last integral can be bounded from below and from above by integrals of the following 



type 



((l + x) 1 -"- 1) +d 



Q 2 exp(-x/2) + (1 + x 
1 + C 3 



,-1 



1 — a 

with some nonnegative constants Cj, i = 1, 2, 3. Integrals in (jlj) have the following asymp- 
totics: 

P = C 3 (j-^ [{d + 1) 1_Q - 1] - C 4 ln(d + 1)^ + 0(1) as d -)• 00 



with some positive constants O3, O4. Therefore, using the asymptotics 

1 



n(tZ) 



1 + a 



^1 + d) a as c? — )■ 00 



(4.10) 



one obtains that the asymptotics of the coarsening law for f)4.9p coincides up to mul- 
tiplicative constants with ( 14. 3 aft already obtained in the example a) for the monotone 
distributions. Consequently, we conclude that as in the example a) the coarsening rates 
are algebraic with power a /{a — 1) for a < 1 and exponential for a > 1. 

A more simple but rather formal proof of this fact is as follows. Let us fix a large 
number A such that asymptotics (14. lUp holds for all d > A with a good precision. The one 
has 



BT(d) 





n(x) 




f d 


n(x) 


/ n(x) In 
/o 


dx + 


/ n(x) In 

J A 


n(d) 


n(d) 



dx 



f d 


n(x) 


/ n(x) In 


n(d) 


J A 



dx. 



The last integral in view of ( I4.10P is of the type considered already in example a). Therefore, 
the term 0(1) + a\n(d + 1) x 0(1) produces no change in the asymptotics of T(d) and, 
consequently, the coarsening law coincides up to multiplicative constants with (j4.3aj) . 
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t x 10* 1 x10* 

Figure 4: Comparison of droplet position evolution obtained from the ODE model (dots) 
and the lubrication model (solid line) in the strong-slip case with e = 0.025, Re = 0, 
/5 = 10. Initial profile of four droplets presented on the top-left plot of Fig. |5]was used. 
Subsequent collisions of the second (left) and the third (right) droplets with the first one 
are shown. 



5 Numerics 



In the last two sections starting from the system fll.2al) - fll.2b[) we first derived a closed 



ODE model fl2.32p - fl2.33p describing coarsening dynamics in an array of initially N + 1 
metastable droplets. Then we looked at its limiting case j3 — > oo described by (12.341) - 
(I2.35P and more precisely on its leading order version as e — > given by (13.11) . Next, we 
found out that for a special initial data satisfying (13.21) one can obtain the explicit solution 
to (13. ip given by (13.41) . Assuming additionally (I3.6P and that the distances in the array are 
ordered decreasingly we derived an explicit coarsening law (13.101) . Finally, we obtained its 
continuous counterpart (I3.14p . In this section we systematically compare numerically the 
solutions of subsequent models in the derived model hierarchy and check coarsening laws 



5.1 Comparison between models 

Here we compare solutions to the full ODE system fl2.32l) - fl2.33l) with those of the strong- 
slip system fll.2ap - fll.2bp and its limiting cases (jl.8p and fll.7ap - fll.7bp as (3 — > and 
(3 — > +oo, respectively. For the solution of PDE systems we used a fully implicit fi- 
nite difference scheme derived and applied already to fll.2al) - (11.2bp and its limiting cases 



in 



15l. 12814301] . The numerical solutions for (I2.32p - fl2.33p were obtained applying a fourth- 
order adaptive time step Runge-Kutta method and using updating rules (I2.38p - fl2.40p after 
each subsequent coarsening event. In the case of PDE system the corresponding pressure 
evolution was calculated using finite-difference discretization of the term H e (h) — d xx h. 

In Fig. H] starting from an array of four droplets we compare evolution of positions 
resulting from PDE and ODE models for two subsequent collisions. Fig. H] shows that 
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Figure 5: Comparison of droplet evolution obtained from the ODE model (dots) and the 
lubrication model (solid line) in the strong-slip case with e = 0.025, Re = 0, and different 
P. Upper row: Initial profile of four droplets (left) and their pressure evolution in the 
intermediate-slip case until collapse of the second droplet (right). Lower row: Starting 
from the same initial profile position evolution until collision of the second droplet with 
the first one for (3 = 5 (left) and (3 = oo (right) is shown. 



the absolute deviation between results stays uniformly 0(e) also after subsequent collision 
events. 

In Fig. [5] starting from the same array of four droplets we compare solutions for 
different slip lengths 0. In the cases (3 = and (3 = oo we compared solutions to (II. 8p and 
(JEZaJ-flEZE) and those of ( l2T36|) - (l2T37j) and (rOi|) -( l2T35|) . respectively. Again for all f3 the 
absolute deviation between PDE and ODE results is 0(e). Note, Fig. demonstrates the 
fact pointed already in 2l|, |27| that in the intermediate-slip case the coarsening dynamics 
is governed mostly by collapse mechanism while in the strong-slip case with moderate and 
large (3 by collisions. 

Finally, Fig. [6] shows that for arrays being taken initially with sufficiently small pres- 
sures the migration subsystem ( 13. ip approximates with a high accuracy (at least 0(e)) 
the full ODE system ( I2.34p -( 12.35p describing the case (3 = oo. Note, that the migration 
path of the fourth droplet Fig. E]is considerably large ~ 50. Hence, the reduction from 
( I2.34p -( l2.35p to (13. ip does not constrains the initial distances to be small. 
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Figure 6: Comparison of droplet position evolution obtained from the ODE model (|2.34p - 
(I2.35P (dots)with e = 0.025 and its leading order subsystem f l3.ll) (solid line). Left plot: 
The initial profile of five droplets. Pressure of the first four and the last droplet are 0.01 
and 0.001, respectively. Right plot: Collision of the fourth droplet with the largest last 
one. 

5.2 Coarsening rates 

Here using the explicit solution ( 13. 4 p for the system (13.1 p we check numerically the discrete 
and continuous coarsening laws (I3.10p and (13.141) . In Fig. [7] we take k = 20 families of 
initial distances as prescribed in (13. 7p with the corresponding pressures satisfying (13. 2p . 
(13. 6p and order them non increasingly in space. Next, we compare the subsequent times 
for each m-th family (1 < m < k) to be absorbed given on one hand by the analytical 
law (13.101) and by iterative calculation using (I3.4p .( 13"3|) between each subsequent collision 
on the other. Naturally, one finds out the exact coincidence between them. Note, as each 
collision is comprised of an absorption by the largest droplet of a smaller one one do not 
need to update the position and pressure of the former one. This is because its position is 
fixed due to (12.311) to x = L and the pressure to the leading order does not change due to 

flHSD. 

In Fig. [S] we show numerical results for continuous coarsening rates for three initial 
distributions taken from the family (14.11) with different a and one Gaussian distribution 
considered in examples a) and c) of the previous section, respectively. To obtain the 
coarsening rates numerically we first sampled N ^> 1 distances according to the given 
initial distribution. After ordering them non increasingly in the initial configuration we 
substitute them as initial data into (13. ip and solve the latter one iteratively using (13. 4p . 
Note, that due to an extremal simplicity of (13. 4p one can effectively model numerically a 
huge number of droplets N w 10 7 just using capabilities of a personal computer. 

Fig. [8] shows that thus obtained numerical coarsening rates coincide for large times very 
well with the analytical ones prescribed by the law (I3.14p and found out in examples a) and 
c) of the section 4. In the case of (14. ip with a = 1 one has the exact coarsening law (I4.4p and 
therefore a good coincidence for all times. In the case of ( 14. ip with a^lwe compared our 
numerical results with the asymptotic law (14. 5p . while for the Gaussian initial distribution 
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Figure 7: Comparison of the subsequent collision/absorption times for each m-th family 
(1 < m < 20) given by the discrete coarsening law (13.101) (dots) and iterative calculation 
using the solution (13.41) (solid line). Left plot: A part of a typical initial profile under 
consideration. Right plot: plot of the absorption times versus the family number. 

with (14. 7|) . Note that a certain deviation between numerics and analytical laws starting 
at the very end of the considered time interval is caused by a numerical error increase in 
sampling of large distances according to initial probability distributions. 



6 Conclusions and discussion 

In this paper we started from the lubrication equations (ll.2al) - fll.2bj) describing dewetting 
process in nanometric polymer film interacting on a hydrophobically coated solid substrate 
in the presence of large slippage at the liquid/solid interface. This model describes a dis- 
tinguished and important regime within a lubrication scaling. In particular, it incorporates 
as a limiting case of the infinite slip length the well-known model of free films (ll.7ap - (11.7bj) 



studied intensively in applications [la, |3j|. Note, that a similar to (ll.7ap - (11.7bj) system 



appears in the study of viscoelastic threads for which coarsening dynamics of interacting 



droplets was observed also at the experiments, see e.g. [32 . 

Motivated by this we derived the reduced ODE models (I2.32I) - (I2.37I) describing coars- 
ening dynamics of droplets governed by (11.2ap - (11.2bl) and its limiting cases. In the limiting 
case (3 = oo we observed that the migration subsystem (13. ip can be solved explicitly for 
special initial data satisfying (I3.2p . (I3.6p . By (I3.4p - (l3.5p the dynamics of droplets consists of 
sequential collisions of smaller ones with the largest one while their distance is distributed 
uniformly between the remaining drops. 

Similar models were suggested for collapse/collision dynamics of breath figures by Der- 
rida et al. [33[ basing on heuristic arguments. There authors considered the 'cut-in-two' 
and 'past-all' models where the distance of the smallest droplet was divided between two 
neighbors or pasted as a whole to one of them, respectively. The breath figures of [33[ found 
later interesting analogs in the reduced coarsening models arising from the Allen-Calm and 
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Figure 8: Comparison of numerical coarsening rates using sampling of the initial data with 
N 1 and subsequent iterative calculation using (13. 4p (solid line) with those provided 
by the analytical law (13.14p (dots). Upper row: loglog and semilog plots for the initial 
distribution (14. ip with a = 1/2 (left) and a = 1 (right), respectively. Lower row: semilog 
plots for (14. ip with a = 20 (left) and for the Gaussian initial distribution (right). According 
to (14. 5 p and (14. 7p in the chosen axe scales the dots reproduce the linear functions except 
for the upper-right plot where they represent the law (I4.4p 
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Ginsburg-Landau equations, see e.g. a review article 34|. A recent generalization of these 



models and rigorous analysis of their self-similar solutions can be found in [35| . Note that 



(13. ip can be classified due to (I3.5P as a 'cut uniformly between all' type model. We are not 
aware if any heuristic or rigorous analog of it was considered so far in the literature. 

Remarkably, our model (13. ip appears as a reduction of a complicated dynamics gov- 
erned by a high-order lubrication system (jl.7ap - (11.7bl) . Moreover, in contrast to stochastic 
in nature models of 33] if the initial distances in ( 13. ip are ordered non increasingly then 
they coarse in a deterministic and exactly solvable scenario. For the latter case we derived 
the coarsening laws (I3.10P and (13.141) analytically and confirmed them numerically. Inter- 
estingly, the derived law ( I3.14p have a form similar to the Shannon's entropy with respect 
to a certain normalization of the droplet number function n(x). The explanation of this 
fact from the statistical point of view will be presented elsewhere. 

Surprisingly, in contrast to the coarsening dynamics governed by reduced ODE models 
arising from ( II. ip which obey always the law ( II. lip our simple model (13.11) can repro- 
duce any algebraic coarsening rates between zero and infinity as well as exponential ones. 
Moreover, for a family of initial distributions (14. ip we showed existence of a threshold for 
their decay at infinity at which the corresponding coarsening rates switch from algebraic 
to exponential ones. Note, that a similar situation was accounted recentl y fo r self-similar 



solutions to Smoluchowski coagulation equation with certain kernels, see [36|, |37 . 

In view of the above observations it would be natural to extend our deterministic 
collision/absorption model to its stochastic variant withdrawing the non increasing order 
of the distances and thus allowing collisions of random droplets with the largest one. As 
in [33| one could probably look for self-similar solutions of the mean-field approximations 
for thus arising stochastic collision models. A further generalization of the model could be 
a withdrawal of the constraints on the initial data (13. 2p . (13. 6p and thus allowing droplets to 
collide and collapse inside of the domain. Note, that an additional difficulty to handle the 
pressure and position update according to the coarsening rules (I2.38p - (l2.40p would appear 
then. 

Finally, it could be possible to derive two-dimensional analogs of the reduced ODE 
models (I2.32I) - (I2.37I) describing physically coarsening of three dimensional droplets on a 
plane substrate. The two dimensional reduced ODE models arising from (11. ip were derived 
in 20|, |2l| . In 20( a mean- filed approximation for the fluxes between droplets was suggested 
under an assumption of well separation of droplets that is unfortunately not suitable for the 
modelling of droplet collisions because the distance between them tends to zero then. In 
this case one should face a problem of solving a Laplace equation (counterpart to equation 
(12. lip in two dimensions) in a complex domain between droplets occupied by the PL 
region. We expect the same problem to appear for the two-dimensional reduced ODE 
models corresponding to (II .2aj) (11 . 2b[) . 
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A Integral I 

Here we show that integral I defined in (I2.16P converges and integrate it explicitly. Chang- 
ing variable in (I2.16P according to the explicit solution (12.5aj) in the CL region, and using 
matching conditions ( 12. 4 p one obtains 



-oo H\ z ( Hi ) ji 



U'(H) 2(U(H) - U{1)) 
H 2 H 3 



dH 



-5/3 + 2H - 1/3H 3 ,„ 

— . = dH . 

y/2/3- H + H 3 /3H 9 / 2 

Let us make a further change of variables t = 1/H and integrate / explicitly as follows. 

1 -5/3t 4 + 2t 3 - 1/3 , f 1 5t 3 - t 2 - t , 1 
dt= dt. - 



o ^2/3t 3 - 1 2 + 1/3 Jo VQtT3 35(3+^)' 

B Connection between discrete and continuous 
coarsening laws 

Here we show that the discrete coarsening law fl3.10p can be recovered back from (13.141) if 
the initial distribution f(x) has the form (I3.15P i.e. if it is represented by k e N families 
as in ( 13. 7p where we denote 



In this case (I3.13P implies 



n(d) = 1 - } j i' m if de [d m , d m -x) 
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Substituting the last expression in (I3.14p implies 

1 



BT(d T , 



4 In 



+ (dm+i ~ dm) In 



1 _ V i' 

i-E 



+ (4-i - 4) In 



l-i'u 



p=m+l m 



1 _ v fc V 

/—ip=m m 



1 - V %' 
Z-^p=m m 

BT{h m+l ) + 



+ ... + 



+ [Nd m + J2( d p- d 



m jip I In 



p=m 



i-E 



k %' ' 

p=m+l in 



1 - V V 

/sp=m rr, 



Therefore, one obtains recursively that 



p=m 



p-p 



i - v fc a 

1 - V fc 2 V 

Z-^p=m m 



(B.1) 



On the other hand dividing (I3.10p by and proceeding to the limit iV — > oo with /c fixed 
one obtains exactly (IB. 1|) . 
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